A Two-Dimensional Effective Model Describing Fluid-Structure Interaction in Blood Flow: Analysis, Simulation and Experimental Validation

We derive a closed system of effective equations describing a time-dependent flow of a viscous incompressible Newtonian fluid through a long and narrow elastic tube. The 3D axially symmetric incompressible Navier-Stokes equations are used to model the flow.Two models are used to describe the tube wall: the Navier equations, for the curved, linearly elastic membrane and the curved, linearly elastic Koiter shell model.We study the behavior of the coupled fluidstructure interaction problem in the limit when the ratio between the radius and the length of the tube, ε, tends to zero. We obtain the reduced equations that are of Biot type with memory. An interesting feature of the reduced equations is that the memory term explicitly captures the viscoelastic nature of the coupled problem. Our model provides significant improvement over the standard 1D approximations of the fluid-structure interaction problem, all of which assume an ad hoc closure assumption for the velocity profile. We performed experimental validation of the reduced model using a mock circulatory flow loop assembled at the Cardiovascular Research Laboratory at the Texas Heart Institute. Experimental results show excellent agreement with the numerically calculated solution. Major application includes blood flow through large human arteries. Résumé Nous obtenons un système fermé d’équations efficaces, décrivant l’écoulement non-stationnaire d’un fluide newtonien incompressible visqueux à travers un tuyau élastique long et de faible épaisseur. Pour modeler l’écoulement, nous utilisons le système de Navier-Stokes 3D axisymétrique et incompressible. Deux modèles sont employés pour décrire la paroi latérale: les équations de Navier pour une membrane courbe élastique linéaire, et ensuite le modèle de Koiter, d’une coque courbe, élastique linéaire. Nous étudions le comportement du système lorsque le rapport ε, entre l’épaisseur caractéristique et la longueur du tube, tend vers zéro. Nous obtenons les équations efficaces, essentiellement 1D, qui sont du type de Biot avec de la mémoire. Une caractéristique intéressante des équations efficaces est que le terme de mémoire capture explicitement la nature viscoélastique du problème couplé. Notre modèle efficace fournit une amélioration significative par rapport aux modèles 1D standards de l’interaction fluide-structure, qui nécessitent une formule de fermeture pour la vitesse, proposée ad hoc. Nous avons effectué la validation ∗Department of Mathematics, University of Houston, 4800 Calhoun Rd., Houston TX 77204-3476, United States ([email protected]). This research was supported by the Texas Higher Education Board, ARP grant #003652-01122001 and by the National Science Foundation under grants DMS9970310 and DMS0245513 †Corresponding author, fax: 0437287415 , tel: 0437287412 ‡Institut Camille Jordan, UFR Mathématiques, Site de Gerland, Bât. A, Université Claude Bernard Lyon 1, 50, avenue Tony Garnier, 69366 Lyon Cedex 07, France ([email protected]). §Department of Mathematics, University of Zagreb, Bijenička 30, 10000 Zagreb, Croatia ([email protected]). expérimentale du modèle réduit en utilisant la bouclé d’écoulement simulé au Cardiovascular Research Laboratory, Texas Heart Institute. Les résultats expérimentaux montrent un accord excellent avec la solution calculée numériquement. L’application principale inclut l’écoulement sanguin à travers de grandes artères du corps humain.